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Abstract. We have developed a new technique for weak 
lensing analysis, with which the effect of the point spread 
function (PSF) on small galaxy images can be corrected 
for accurately. Rather than relying on weighted second 
moments of detected images, which we show can leave 
residuals at the level of a percent in the shear, we di- 
rectly fit (stacked or individual) galaxy images as PSF- 
convolved, sheared circular sources. We show by means of 
simulations that this technique is able to recover shears 
well below the percent level for a variety of PSF shapes, 
and that its noise properties are similar to existing meth- 
ods. 

Key words: Methods: data analysis - gravitational lens- 
ing 



1. Introduction 

Gravitational lensing is one of the most powerful and di- 
rect methods for studying the gravitational potentials of 
massive objects in the universe. An important type of lens- 
ing study is 'weak lensing.' It is the study of mild system- 
atic distortions of background sources as their light rays 
are perturbed by gravitational fields on their way to us. 
Weak lensing has already provided importan t resu lts in 
the study of ga laxy clusters (e. g., Ty son et al. 1990; Bon- 



net et al. 1993; Fahlman et al. 1994; Squires et al 



Fischer et al. 1997; Clowe et al. 1998; Hoekstra et al 



J1997; 
19981) , 



halos of individual galaxies (e.g., Brain erd et al. [1996| ), and 
large-scale structure (Schneider et al. 1998). As the tech- 
niques are becoming better understood, research is pro- 
gressing to the search for weaker and weaker distortions, 
which would enable the outer regions of galaxy clusters 
and galaxy halos, as well as lensing signals from large- 
scale structure (e.g., Jain & Seljak 1997 ; Kaiser 1998|) , to 
be studied. 

To be able to detect such very weak signals, it is im- 
portant to accurately remove the dominant systematic ef- 
fect affecting weak lensing measurements: anisotropy of, 
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and smearing by, the point-spread function (PSF). In 
this paper we will first investigate the limits of the most 
commonly-used technique for weak lensing analysis, de- 
vised by Kaiser et al. ( 1994 henceforth KSB). We will 
show that after PSF anisotropy correction, residual effects 
on the order of 1% shear are difficult to avoid with this 
method, even for moderately elongated PSF's. Since the 
ability to detect percent signals is important for a variety 
of scientific questions, we have therefore devised a new 
method which does not have such residuals, but which 
nevertheless has noise properties comparable to those of 
the KSB method. 

There are several other methods for PSF anisotropy 
correction in the literature. The Autocorrelation Function 



method of Van Waerbeke et al. ( 1997 ) is a variant of the 
KSB method in which not individual galaxy images, but 
the autocorrelation function of many of them, is analyzed. 
The Bonnet & Mellier ( 1995 ) method uses a different aper- 
ture weighting function from KSB, and treat s the PSF 
convolution as a shear term. Fisher & Tyson (1997) con- 
volve the image with a kernel constructed to make the 
PSF rounder again. A more sophisticated such kernel has 
recently been presented by Kaiser (1999). 



2. Limitations of the KSB formalism for PSF 
anisotropy correction 



We first examine the Kaiser et al. (1995) method, following 
our earlier limited investigation in the context of analysis 



of Hubble Space Telescope images (Hoekstra et al. 1998, 
Appendix D). 

2.1. The KSB method 

The technique of weak (or statistical) lensing involves 
measuring the systematic, gravitationally induced, distor- 
tion of background images behind a gravitational lens. In 
the weak lensing regime small background images are dis- 
torted by a shear (71, 72) and a convergence k, whose com- 
bined effect is represented by the mapping 
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Fig. 1. An example of a PSF with radially varying ellip- 
ticity, but zero overall polarization. (Eq. ^ with 6 — 0.3). 
Contours differ by a factor of 2^/^. The dotted curves are 
circles, shown for comparison. 



where gi = 7i/(l — k). For simplicity, in what follows we 
neglect k as it is small in the weak lensing regime, and 
pretend we are deriving the true shear 7 instead of the 
reduced shear g. Thus, our results on shape measurements 
are valid, but their interpretation as a lensing signal may 
require consideration of the (1 — k) factor. Our analysis 
makes no assumptions on the smallness of gi, though. 

Kaiser, Squires and Broadhurst (1995, henceforth 
KSB) describe a method for recovering (71,72) from im- 
ages of distant galaxies. Essentially, they derive galaxy el- 
lipticities from weighted second moments of the observed 
images, and then correct these for the effects of the weight 
function and of smearing by the point spread function 
(PSF). By averaging over many galaxies, which are as- 
sumed to be intrinsically randomly oriented, the effect of 
individual galaxy ellipticities should average out, leaving 
the systematic lensing signal. The KSB method has proved 
to be very effective, especially in the study of galaxy clus- 
ter potentials. 

KSB define various "polarizabilities" , which express 
the ratio between an input distortion (gravitational shear 
or PSF anisotropy) and the measured polarization 

\ Jw I Tw ' Tw I Tw J ^ ^ 
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Fig. 2. The polarization of the PSF of Figure yj as a func- 
tion of the weight function's Gaussian radius w. The plot 
shows that with compact weight functions such as those 
that are required to suppress photon noise, even mildly 
elongated PSF's may erroneously yield a polarization of 
several percent. 

of an image f{x,y). These polarizations are derived from 
weighted second moments 

= J f{x,y)x^Wir)dxdy, etc., (3) 

of the image intensities, where W is a weight function 
which goes to zero at large radii. The weight function 
is required as otherwise the sky noise in the outer parts 
of the image dominates the measured moment. The sig- 
nificance of the measurement is optimized by taking the 
weight function to be relatively compact, of a size compa- 
rable to the image itself. 

Details of the method can be found in KSB, and 
in Hoekstra et al. (1998, henceforth HFKS), where a 
few small errors in the formulae of KSB were corrected. 
For the purposes of the present paper, it is sufficient to 
know that in the KSB formalism, the "smear polarizabil- 
ity" P^™ defines the ratio between the PSF anisotropy 
P = [Ixx — Iyy,'2Ixy), Constructed from the unweighted 
second moments lij of the (normalized) PSF, and the re- 
sulting change in image polarization e. The "shear polariz- 
ability" P^^ is the ratio between the applied shear (71, 72) 
and the resulting change in the image polarization. KSB 
show how the polarizabilities can be derived from higher 
weighted moments of the observed PSF and galaxy im- 
ages. 
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2.2. How accurate is KSB? 

In the context of ground-based cluster weak lensing, the 
KSB method works well. Nevertheless, it does involve 
some approximations. Now that weaker and weaker sig- 
nals are of interest, it is therefore important to understand 
the limitations of the method. As already discussed by 
HFKS, for strongly non-Gaussian PSF's the KSB method 
does not completely correct PSF anisotropy. This is par- 
ticularly true when analyzing small galaxies in deep HST 
images, where it turns out that the choice of weight func- 
tion in eq. ^ is important. 

A simple PSF model can be used to illustrate why 
such residuals are, at some level, unavoidable. Consider 
the following PSF: 



(4) 



where G{a^,b^) is a unit-area Gaussian of x- and y- 
dispersions a and b. d is a small parameter. The case 
i5 = 0.3 is plotted in Figure |l]. 

The PSF of equation || has exactly zero anisotropy p: 
the second moments in x and in y are equal. However, 
the ellipticity of the PSF varies with radius, which means 
that the weighted second moments are not equal: weighting 
the central parts more will enhance the a;-moment prefer- 
entially. In fact, it is easy to show that the polarization 
constructed with weighted moments is 0{S). The precise 
result for a Gaussian weight function W = cxp(— ir^/w^) 
is 



ei 



9 {7 + 5w'^ + w^) S 



2 (1 + u;2) (4 + ^„2) (20 + 16 ^2 + 5 



The polarization of the 6 ~ 0.3 PSF plotted in Figure Q 
is plotted in Figure ||. It has a value of 0.03 near w = 2, 
roughly the radius of maximum significance which should 
be used to minimize photon noise in the polarization mea- 
surement. 

The KSB polarizabilities are derived assuming that 
the PSF can be written as the convolution of a com- 
pact anisotropic part with an extended circular part (KSB 
eq. Al). This assumption allows the anisotropy to be char- 
acterized in terms of p only. However, our example shows 
that this assumption may be too restrictive: it effectively 
couples the radial intensity profile of the PSF with its el- 
lipticity profile. For example, a single Gaussian with con- 
stant ellipticity can be written as such a convolution, but 
a sum of two elliptical Gaussian such as the PSF of eq. || 
cannot. The systematic errors that arise are the result of 
this. 

3. A new method 

Here we present a new method, with which the PSF effects 
can be corrected for with greater accuracy. The essence 
of the method is not to work with the moments of the 
observed images; instead each image is fit directly as a 



PSF- convolved, sheared circular source of unknown radial 
profile. 

Assume for the moment that we have managed to sum 
the images of many galaxies into an 'average galaxy' im- 
age g{x,y). Analysing a stacked galaxy image is similar 
to the approach discussed by Lombardi & Bertin ( 1998 ), 
who average image second moments before corrections are 
applied. It differs from methods such as KSB or Bonnet & 



Mellier (1995) in which galaxies are individually corrected 
for PSF effects before they are combined to produce a 
shear estimate. 

Intrinsically, 'g is circular if the galaxies are randomly 
oriented, but the image we observe has been distorted first 
by gravitational lensing shear, then by the atmospheric 
seeing, and finally by the camera optics. The observed g 
is therefore a sheared circular source, convolved with a 
(known) PSF. We therefore fit g directly to such a model, 
with the minimum of further assumptions. This approach 
addresses the apparent difficulty in the KSB methodol- 
ogy in the case of radially changing ellipticity profiles: a 
sheared circular source has constant ellipticity at all radii, 
and so after convolution with the PSF only a subset of 
"allowed" ellipticity profiles remain. 

Assuming that the PSF is known, e.g., from analysis 
of star images in the field, the model for 'g is specified by 
an unknown radial brightness profile, and by the shear 
parameters (71,72) that we are interested in. In prac- 
tice we model the radial profile as the superposition of 
several Gaussians of different fixed widths, and unknown 
amplitude. We have found that the following recipe for 
assigning the basis functions gives good results: (i) de- 
termine the best-fit circular Gaussian radii to the ob- 
served PSF and galaxy images, rpsF and tqal- (ii) Take 



,2 

GAL 



' PSF 



1^/2 as an estimate for the intrinsic ra- 



dius of g. (iii) Use four components to describe the radial 
profile of g, with Gaussian radii (0.5, 1, 2, 4) x r. 

The algorithm is laid out in Figure ||. We now de- 
scribe the results of tests to verify the accuracy of the 
PSF anisotropy correction, and to investigate how well it 
fares in the presence of noise in the images. 

3.1. Simulations in the absence of noise 

We tested how well PSF anisotropy can be corrected for by 
considering the case where there is no gravitational shear, 
only a range of PSF shapes of varying anisotropy. An ac- 
curate analysis should yield zero shear after correction for 
the PSF. On a large number of model images, described 
below, we compared the results of the algorithm of Fig- 
ure 1^ with those from the KSB algorithm as described in 
HFKS (implying in particular that the same weight func- 
tion is used in the derivation of polarizations and polariz- 
abilities of galaxy and PSF images). The weight function 
was taken to be the best-fit circular Gaussian to the post- 
seeing galaxy image. 
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Fig. 3. The schematic algorithm used to derive the shear from an observed mean galaxy and PSF image. 



The algorithm described in this paper directly yields 
an estimate for the shear. In the comparisons, the KSB 
galaxy polarization after seeing anisotropy correction was 
divided by the "pre-seeing shear polarizability" , for 
which we use the expression given by Luppino & Kaiser 
(1997). 



3.1.1. Double-Gaussian images and PSF 

In most of our simulations, we modeled the round average 
galaxy images as 



g^G{rg) + G{kgrg) 



(6) 



where G{a) is a unit-integral Gaussian of dispersion a (a 
double-Gaussian PSF was also considered by ???REF???). 
The parameter kg is unity for a Gaussian profile, and is 
larger for more radially extended profiles. A reasonable, 
though admittedly crude, approximation to an exponen- 
tial profile is given by setting fcg = 2, while fcg = 3 gives 
a reasonable approximation to a de vaucouleurs profile 
(Figure ||). 




Fig. 4. The double- Gaussian profiles used in the modeling 
in this paper. Left, the k = 2 profile is plotted logarithmi- 
cally to show its similarity to an exponential profile; right 
the fc = 3 profile is plotted logarithmically vs. r^^^ to show 
it is similar to a de Vaucouleurs model. 

The PSF's were modeled in a similar way, but with 
anisotropy. Writing now G{a, b) for a Gaussian with x- 
and y-dispersions a and b, we have 

PSF = G (rp, (1 - ei)rp) + G {kpTp, (1 - e2)kprp) . (7) 
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Again we introduced a shape parameter kp, but we also 
included ellipticities for the two components. We con- 
sidered three kinds of PSF ellipticity profile: we either set 
ei — 62 (constant ellipticity with radius), or set one of the 
e's to zero, to give a radial increase or decrease of the PSF 
anisotropy. These three possibilities, though by no means 
exhaustive, form a representative set of PSF's. 

The advantage of the multiple-Gaussian formulation 
is that the PSF convolution can be done analytically. We 
thus constructed a large number of PSF-smeared galaxy 
images, calculated the polarizations and polarizabilities 
following HFKS, subtracted the PSF anisotropy correction 

to the galaxy polarization, and divided the result by Lup- 
pino and Kaiser (1997) 's pre-seeing shear polarizability 
P^. We then compared these results with the results of 
our implementation of the new fitting algorithm of Fig- 
ure ||. 

The results of the simulations are presented in figures 
|, I and §. They show that the KSB method can suffer 
from systematic residuals around the 0.01 shear level once 
the PSF ellipticity exceeds 0.2 or so, whereas this is not 
so for the new method developed here. The systematic 
effects are strongest for small galaxies, for PSF profiles 
with long tails, and for radially increasing PSF ellipticity. 
The effect is clearly driven by the PSF shape, not by the 
galaxy brightness profile. 

Notice that in the constant-ellipticity case (Figure 
with a Gaussian PSF (fcpsf = 1) the residuals left by the 
KSB method are high order in PSF ellipticity, but that for 
non-Gaussian PSF's a low-order residual dominates. (We 
have verified this result analytically using symbolic math- 
ematics.) This is a consequence of the fact that only the 
single elliptical Gaussian PSF can be written as a convo- 
lution of a compact anisotropic function with a round ex- 
tended one, as assumed in the KSB derivation. It is clearly 
important to test algorithms not only for single-Gaussian 
PSF's! 



3.1.2. A WFPC-2 PSF 

In order to test whether our results are specific to the 
double-Gaussian formulation of the PSF, a test was also 
performed with a model PSF for the WFPC-2 camera on 
the Hubble Space Telescope. The model was generated 
with the Tiny TIM software package, provided on-line at 
STScI by J. Krist. An oversampled PSF was calculated 
for a position near the corner of CCD#4, and convolved 
with a Gaussian circular galaxy of FWHM 0.25arcsec. 
This 'galaxy' and the PSF (Figure ||) were then binned 
to a resolution of half a WFPC-2 pixel to avoid under- 
resolving the PSF, and analyzed as above. The results are 
summarized in table ^ and confirm the results obtained 
from the large number of double-Gaussian simulations de- 
scribed earlier. 
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Fig. 8. PSF and simulated galaxy image for an observa- 
tion in the corner of one of the WFPC2 CCD's. Axis units 
are 0.05 arcsec. 

Table 1. Results of a simulation based on a WFPC-2 
PSF calculated using the TinyTim software. In agreement 
with earlier results (HFKS), the KSB technique appears 
to over-correct for the anisotropic WFPC-2 PSF images 
slightly. 



Weight function radius: 
PSF polarization: 
Uncorrected galaxy e: 
PSF-corrected galaxy 7: 



1.3 WFPC-2 pixels 
(-0.078, 0.024) 
(-0.020, -0.006) 
(-0.000, -0.019) (KSB) 
(-0.000, -0.001) (new method) 



3.2. Noise properties 
3.2.1. Analytic estimate 

The error on the estimated shear due to photon noise can 
be estimated as follows. Let the 1-a error on each pixel of 
g he a (for simplicity we take this to be the same on every 
pixel, appropriate for background- limited work) . Then the 
fit involves finding the minimum of 



X 



(8) 



where P{x,y) is the PSF, /(r^) is the intrinsic radial pro- 
file of the average galaxy, (g) denotes convolution, is the 
position of the fcth pixel, and T is the distortion matrix 
of equation ^ If the fit parameters 7^ are uncorrelated 
with the radial profile, their inverse variances are given by 
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Fig. 5. The result of correcting simulated unsheared images for PSF anisotropy, following the KSB method (solid lines) 
and the method presented here (dashed lines). S'yi is in each case the shear that is deduced after the PSF correction, 
and should be zero for a perfect analysis. The fc's are luminosity profile shape parameters for galaxy and PSF, and 
are explained in the text. In each panel the heavy line represents the case where the galaxy image is intrinsically of 
the same radius as the PSF (r^ = rp in eqs. ^ and 0), the lighter lines those where the galaxy is 0.5 (upper) and 
1.5 (lower) times this size. The PSF ellipticity is constant with radius in these simulations. While for the fcpsp = 1 
case (Gaussian PSF) the KSB method leaves a residual which is third-order in PSF ellipticity, other PSF luminosity 
profiles give rise to first-order residuals. The residuals of the new method in the upper panels disappear if more radial 
components are used in the fit for highlighting that the dominant source of error in this method is the extent to 
which the radial profile is modeled correctly. 
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Fig. 6. As figure ^, but only the outer component of the PSF is eUipticaL 

\d-^X^ /d^i- For example, at the best fit and respectively. Then the l-cr error on 71 evaluates to 



where we have used the results that the PSF-fitting error 
The right-hand side of equation 1^ can be estimated assum- on F for a Gaussian source is 5F = 27ri/2r.GALO-. (The 
ing that the PSF and observed average galaxy are Gaus- error on 72 is the same.) We have verified this formula 
sians with dispersions rpsF and tqal pixels of integral 1 by means of simulations, similar to those described below. 
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Fig. 7. As figure ^, but only the inner component of the PSF is elliptical. 



Equation 10 shows the expected increase in noise for small 
objects, as well as the limit, even for large objects, of 



o'(7i) ^ 



6F_ 



(11) 



3.2.2. Simulations 



We have checked the sensitivity to noise in the images 
by Monte Carlo simulation. Many realizations of random 
Gaussian noise superimposed on a PSF-smeared, intrinsi- 



cally round galaxy image were analyzed with both algo- 
rithms, and the distributions of the resulting (71,72) es- 
timates compared. Selected results are shown in Table ^ 
Interestingly, the dispersions in the shears derived with 
both methods are very similar over a range of galaxy sizes. 
As we have already seen, the small bias in the results from 
KSB is present in the simulations with non-circular PSF's, 
but not in the method advocated in this paper. 

A possible way to avoid the systematic residuals of 
the KSB method is to increase the radius of the weight 
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function W in eq. ^, since the problems arise from the 
imperfect way in which the polarizabilities represent the 
effect of W. However, the primary function of W is to 
control the noise in the images. Doubling the Gaussian 
radius of W does in fact improve the anisotropy correction 
in the mean, but at the cost of almost doubling the noise 
on the result (see Table ^) . 

A by-product of our algorithm is an estimate of the 
intrinsic radial profile of g. In practice, this estimate ap- 
pears to be rather sensitive to the noise, especially for 
small images — not surprising given that this is effectively 
a deconvolution, albeit a constrained one. Nevertheless, 
it may be possible to use the information in the best-fit 
radial profile in several ways. If a suitable prior for the 
intrinsic radial profile of the average galaxy selected can 
be formulated (e.g., by combining results over a wide field, 
or from deeper, higher resolution images), this information 
might help to refine the best-fit shear solution further. Al- 
ternatively, the width of g might be used to attempt to 
derive the lensing convergence k directly, since in principle 
it is a direct measure of the magnification of faint galax- 
ies. This possibility is yet to be explored in detail, but is 
likely to be difhcult in practice. 

3.2.3. The effect of centroiding errors 

The centroid of an image can be determined in differ- 
ent ways, each of them susceptible to errors due to pho- 
ton noise. The effect of centroiding errors on the summed 
galaxy image will be a convolution with the distribution of 
centroid errors. Thus, the PSF needs to be convolved with 
this distribution before analysis of so that the effect of 
the centroiding error can be compensated. 

4. Galsixy-by-galaxy application 

The method as described so far involves analysing the av- 
erage galaxy g. Very accurate shear measurements require 
g to be the average of a large number of galaxies (~ 1000 
for a 1-(T shear accuracy of 0.01), otherwise intrinsic ellip- 
ticity scatter will dominate. However, constructing such a 
g is only possible if the shear and the PSF are constant 
over a large part of the image. Often this is not the case. 

To cope with this limitation, we have therefore exper- 
imented with the algorithm in 'galaxy-by-galaxy' mode, 
where the algorithm is applied to individual galaxies and 
the resulting shear estimates averaged. Mathematically 
this approach is not perfect, because it involves fitting 
a constant-ellipticity model to individual galaxies even 
though this is not necessarily appropriate. Nevertheless it 
turns out to work better than might have been expected, 
and better than existing methods. 

We tested this approach on various model galaxies, 
of differing axis ratios. To simulate typical galaxies, we 
include a round, central 'bulge' component, and an outer 
'disk' of axis ratio between 0.1 and 1. (Simulations with 



different bulge axis ratios yielded the same results.) These 
were placed at all orientations, smeared with an elliptical 
PSF, and analysed with the algorithm described above. 
The best-fit (71,72) values thus derived for each galaxy 
are then averaged to give an estimate of the shear. 

As may be seen in Figure ^, the algorithm performs 
very well, essentially correcting all PSF anisotropy signal 
in the measured shear. By comparison the slightly biased 
answer returned by the KSB algorithm is apparant as be- 
fore. Residual systematics of the new method are at the 
level of a few tenths of a percent. 

5. Summary 

In this paper we have studied possible systematic errors 
arising from the correction for anisotropic point-spread 
functions in weak lensing analyses based on the well- 
known Kaiser et al. (1995) method. While such effects 
are small, generally below a few percent in the deduced 
gravitational shear components (71 ,72), they are at a level 
that is important for studies such as galaxy-galaxy lensing, 
lensing by large-scale structure or cluster lensing at large 
radii. A range of simulations shows that modelling the 
PSF as a convolution of a compact anisotropic function 
with a more extended, circular function, which underlies 
the KSB formulation, is not sufficiently general to describe 
many PSFs, and leads to these systematic residuals. 

We have presented a new algorithm with which to 
carry out the PSF-correction in a single fitting step, and 
show with simulated images that the low-level residuals 
left in the KSB analysis can thus be avoided. The whole 
image is used in the fitting, so that not just the lowest mo- 
ments are used to characterize the image shapes. We have 
also shown that the noise properties of this algorithm com- 
pare well with those of KSB. While mathematically the 
algorithm requires an intrinsically circular source, as may 
be constructed by stacking many observed galaxy images, 
in practice nearly unbiased results can also be obtained 
when the algorithm is used to correct individual galaxy 
shapes for the PSF. This galaxy-by-galaxy application of 
the algorithm allows observations with spatially varying 
PSF and/or shear fields to be handled. 

Weak lensing is a unique technique with which to study 
gravitational potentials at large radii in galaxy clusters, 
galaxy halos and in the field. The present method holds 
the promise of allowing a little more information to be ex- 
tracted from the large volumes of data that will be gath- 
ered with the coming generation of wide-field imagers. 

Acknowledgements. I would like to thank Peter Schneider, 
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ings of the manuscript and for suggesting several improve- 
ments. 
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Fig. 9. The derived shear values from the algorithm applied to individual elongated 'disk-fbulge' galaxy images, after 
smearing with an elongated PSF. Each ring represents a galaxy of different shape, seen at many orientations. Left: 
the raw (61,62) polarizations measured with the standard KSB method, without correction for the PSF. Centre: 
the result of applying the KSB prescription for PSF anisotropy and circularization correction. The small bias seen 
before remains. Right: the result of the new algorithm on the same galaxies. In the latter case, the correct average 
shear estimate (zero) is recovered even though individual galaxies are not correctly described as intrinsically constant- 
ellipticity sources. The three rows refer to the three kinds of PSF ellipticity profile considered in figures |[-^: constant 
with radius (top), outward-increasing (middle), and outward-decreasing (bottom). 



Konrad Kuijken: Weak weak lensing 



11 



Table 2. Results from representative noise simulations of the KSB method and the one presented in this paper. In each 
case, 100 noise realizations (noise per pixel of 0.001, with rpsf = 2pixels, fcgai = 2 and total flux 4) were analyzed with 
the standard KSB method, with the KSB method using a weight function double the radius of the best-fit Gaussian, 
and with the new method described in this paper. The first six simulations were of cases without PSF anisotropy, and 
in the last six the PSF has a constant axis ratio of 0.7. In all cases, the dispersions of the standard KSB method and 
the new one are very similar, but note the imperfect correction from the KSB method. The simulations with a larger 
weight function show that the PSF anisotropy is corrected better, but at the cost of increased noise. 
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